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ABSTRACT 


A dynamic analysis of permanent magnet stepping motors for multistep operation is 
described. Linearized transfer functions for single-step responses are developed. Mul- 
tistep operation, when the motor is driven by a current source, is analyzed by using 
phase plane techniques. Failure of the motor to operate for fixed stepping rates and 
load torques is discussed. Dimensionless curves showing maximum stepping rate as a 
function of motor parameters and load torque are developed and experimentally verified. 
The curves allowed the maximum stepping rate of a motor with viscous, inertial, and 
torque loads to be predicted from simple measurements of the motor parameters. 



DYNAMIC ANALYSIS OF PERMANENT MAGNET STEPPING MOTORS 

by David J. Robinson 
Lewis Research Center 

SUMMARY 

Permanent magnet stepping motors are being applied as digital actuators. Analytical 
models of stepping motors have been developed that analyze the dynamic response to a 
single step. They have not, however, been expanded to analyze stepping motor perform- 
ance during multistep operation. 

An analytical study is presented herein that examines stepping motor performance 
during both single-step and multistep operation. A linearized single-step model is 
developed which allows the stepping motor to be defined in terms of a natural frequency 
and a damping ratio. A nonlinear analysis is also developed which, when a constant cur- 
rent source is assumed, allows the stepping motor to be analyzed for fixed stepping rates 
and applied load torques. It was found that the permanent magnet stepping motor cannot 
respond to a step command when the applied load torque becomes greater than 0. 707 of 
the motor's stall torque. Also the motor cannot follow a sequential set of step commands 
if, during the sequence, the rotor lags the command position by more than two steps. 

Based on these results, we developed and experimentally checked a set of dimension- 
less curves which express maximum normalized stepping rate as a function of normalized 
damping and normalized load torque. These curves allow the maximum stepping rate of a 
motor with viscous, inertial, and torque loads to be predicted from simple measurements 
of the motor parameters. 


INTRODUCTION 

In recent years, the use of digital control has found increased application in the field 
of automatic control systems. One inherent problem in applying digital control is in se- 
lected suitable digital actuators. One promising actuator is the permanent magnet (PM) 
stepping motor. 

The application of stepping motors is not new. Proctor (ref. 1) reports that stepping 
motors were first used in servomechanisms in the early 1930's. During the development 



of closed-loop servomechanisms in the World War n years, proportional analog activa- 
tors largely replaced stepping motors in actuator applications. With the space age and 
the development of automatic digital systems, problems arose in using analog closed- loop 
servomechanisms because of the need for digital- to- analog (D/A) conversion techniques. 
This problem was particularly important in space applications because of the extra weight 
required for the D/A equipment. 

The stepping motor has found increased applications as a digital actuator because it 
does not require D/A conversion. Nicklas (ref. 2) utilized a stepping motor as an actua- 
tor in a spacecraft instrumentation system. Giles and Marcus (ref. 3) reported on the use 
of stepping motors as control drum actuators for the SNAP-8 program. 

To increase the application of the stepping motor, much analytical work has been 
done. Bailey (ref. 4) compares stepping motors with conventional closed-loop positioning 
systems. The stepping motor was modeled, for a single-step response, as a second- 
order linear approximation. O’Donohue (ref. 5) and Kieburtz (ref. 6) have developed 
similar mathematical models. These particular second -order models cannot be extended 
to adequately describe the dynamics of the stepping motor for multistep inputs. Bailey 
analyzed multistep inputs on a limited basis, by assuming that the time between applica- 
tion of the step command was long compared with the settling time of the response to each 
step input. 

The work described in this report was conducted to develop analytical techniques to 
investigate multistep operation. The analysis is restricted to permanent magnet stepping 
motors; however, the techniques can be extended to other forms of stepping motors, 
electrical, pneumatic, or mechanical. Stepping motor operation is analyzed to determine 
an expression for developed torque. A linearized single-step analysis is presented to ex- 
press motor performance in terms of a natural frequency and damping ratio. Phase plane 
techniques are used to extend the analysis to handle multistep commands supplied from a 
constant current source. From the analysis, a dimensionless curve is developed that al- 
lows the maximum stepping rate of a motor with viscous, inertial, and torque loads to be 
predicted from simple measurements of the motor parameters. 

Simple experimental results from two stepping motors are presented. These results 
are used to check the analytical analysis. 


DESCRIPTION OF PERMANENT MAGNET STEPPING MOTOR OPERATION 

The permanent magnet (PM) stepping motor is an incremental device that accepts 
discrete input commands, and responds to these commands by rotating an output shaft in 
equal angular increments, called steps; one step for each input command. The angular 
position of the output shaft is controlled by the initial location and the number of input 
commands received. The angular velocity of the output shaft is controlled by the rate at 
which the input commands are carried out. 
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Simplified Permanent Magnet Stepping Motor 


A simplified PM stepping motor is shown in figure 1. The PM stepping motor con- 
sists of a stator containing two or more phases wound on salient poles and a permanent 
magnet rotor of high permeability. When a stator winding is energized, a magnetic flux 
is set up which interacts with the permanent magnet rotor. The rotor will move in such 
a manner that the magnetic moment of the permanent magnet will aline with the field set 
up by the stator winding current. In referring to figure 1, assume winding 1 is ener- 
gized such that a magnetic flux is set up with a direction from the face of salient pole 1 
and into the face of salient pole 3. The rotor will aline as shown. If winding 2 is ener- 
gized such that a magnetic field is set up with a direction from the face of salient pole 2 
and into the face of salient pole 4, the rotor will turn such that the south magnetic pole of 
the rotor alines with salient pole face 2. If winding 1 is energized in the manner pre- 
viously described, the rotor will move so that the south magnetic pole is again alined 
with salient pole face 1. However, if the current is reversed in winding 1 such that the 
magnetic field direction is from salient pole 3 and into salient pole 1; the rotor will move 
such that the south magnetic pole alines with salient pole 3 instead of salient pole 1. 

Thus, the position of the rotor of a PM stepping motor can be determined by a discrete 
set of stator winding excitations and a discrete set of current reversals in the stator 
windings . 


Salient 
pole 2 



Figure 1. - Typical representation of permanent magnet stepping motor. 
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Synchronous Inductor Motor 


There are several different PM stepping motor configurations available commer- 
cially. The most common types are derived from two- or four-phase ac synchronous 
motors. The synchronous inductor motor is one of these types and has received consid- 
erable attention in the literature for stepping motor applications (refs. 1 and 6 to 9). 

As an aid in developing some general results, the torque developed by the synchro- 
nous inductor motor is examined. Figure 2(a) shows a simplified cross section of a two- 
phase synchronous inductor motor with four salient poles and five rotor teeth. The rotor 
tooth pitch is 72°, while the salient poles are located every 90°. One step corresponds 
to one- fourth of the rotor tooth pitch, or a rotor movement of 18°. One complete revo- 
lution of the rotor corresponds to 20 steps. Figure 2(b) shows an expanded layout of the 
rotor and stator. The expanded layout is used to simplify the graphical representation 
of the developed torque as the rotor moves relative to the stator salient poles. 

Stepping the rotor is accomplished by reversing the direction of the current in one 
phase while holding the other phase constant. For a two-phase, four- salient- pole motor, 
the salient poles can form four possible magnetic pole combinations: 

(1) Salient poles and Pg^ are north magnetic poles, and salient poles P 1B 

and Pgg are south magnetic poles. 

(2) Salient poles P 1B and Pg^ are north magnetic poles, and salient poles P^ 
and Pg B are south magnetic poles. 

(3) Salient poles P^ B and Pg B are north magnetic poles, and salient poles P^ 
and Pg^ are south magnetic poles. 

(4) Salient poles P 1A and Pg B are north magnetic poles, and salient poles P 1B 
and Pg^ are south magnetic poles. 

Figure 3(a) shows the developed torque as a function of rotor position for combination 
(1) . The resultant curve is the summation of the torque curves due to the individual 
phases and, as will be discussed later, is assumed to be sinusodal. The equilibrium 
points of the resultant curve lie midway between the equilibrium points of the individual 
phases. Figure 3(b) shows the developed torque as a function of rotor position for com- 
bination (2) . The equilibrium points of the resultant curve remain midway between the 
equilibrium points of the individual phases. Comparing the resultant curves of figures 
3(a) and (b) show that the equilibrium points have shifted by 1/4 rotor tooth pitch. In like 
manner, figures 3(c) and (d) show the developed torque for combinations (3) and (4). If 
the motor is stepped by sequentially repeating the magnetic pole combinations, the rotor 
equilibrium points will be shifted by 1/4 rotor tooth pitch for each combination applied. 
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(a) Cross section. 



(c) Expanded rotor-stator layout. 

Figure 2. - Simplified synchronous inductor motor. 
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(a) Stator currents Ij and I 2 : salient poles P 1A and P 2A are north 
magnetic poles,* P 1B and P 2B are south magnetic poles. 
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(b) Stator currents -Ij and I 2 : salient poles P 1B and P^ are north 
magnetic poles; P 1A and P 2 g are south magnetic poles. 
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(c) Stator currents -Ij and -I 2 : salient poles P^g and P 2B are north 
magnetic poles; P| A and P 2A are south magnetic poles. 
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(d) Stator currents Ij and -I 2 : salient poles P 1A and P 2 g are north 
magnetic poles; P|g and P 2A are south magnetic poles. 

Figure 3. - Torque as function of rotor position - both phases energized. 
Arrows indicate direction of developed torque about equilibrium points. 



Analytical Development of Stepping Motor Torque 


Torque is produced on the rotor of the PM stepping motor as the result of an inter- 
action between the flux created by the stator windings and the permanent magnet rotor. 

In the PM stepping motor, the stator windings are wound in coils surrounding each salient 
pole. Assume that an energized stator winding produces a parallel magnetic field of den- 
sity B. Also assume that the rotor of the PM stepping motor is a thin bar magnet. The 
torque relations of a thin bar magnet in a parallel magnetic field are shown in figure 4. 



(b) Maximum and minimum torque developed. 

Figure 4. - Torque relations of thin -bar magnet in parallel 
magnetic field. 


The bar magnet tends to turn in a direction such that its magnetic moment vector M will 
line up with the parallel magnetic field. The torque developed Tp> is the cross product 
of the magnetic moment vector M and the magnetic flux density B: 

T D = M x B (1) 


(All symbols are defined in appendix A. ) When figure 4(a) is used, the developed torque 
becomes 
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T d | = |m| |b| cos a 


( 2 ) 


Figure 4(b) shows that the maximum torque is developed when the permanent magnet is 
perpendicular to the parallel field. When the permanent magnet is parallel with the mag- 
netic field, no torque is developed. 

The magnetic moment of the thin bar magnet is proportional to the pole strength m 
and to the length Z of the bar magnet. Thus, 

|M| = K m Zm (3) 

If the parallel magnetic field is produced by an ideal solenoid of n turns, the magnetic 
density is proportional to the number of turns and to the magnitude of the current I(t) 
flowing through the solenoid turns. 


| B | = Kjnl(t) (4) 

Thus, the torque developed on a thin bar magnet in a parallel magnetic field is given by. 

T D =KjK m ZmnI(t) cos a (5) 

or 

T d = K T I(t) cos a (6) 


where 


K t = KmZmnKj 

Since the maximum torque developed occurs when a = 0, 


T 

max 


= K T I(t) 


Thus, the developed torque can be expressed by 


( 7 ) 


T D = 


T max cos a 


(8) 


Equation (8) may not be an exact model when applied to the PM stepping motor, since 
the stator windings are not ideal solenoids and the rotor is not a bar magnet. The mag- 
netic flux set up by the stator windings has a distribution effect across the surface of the 
salient pole face. Similarly, the rotor has distributional effects due to its configuration. 
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▲ Stable equilibrium points 
• Unstable equilibrium points 


Salient pole Salient pole Salient pole Salient pole 
P 2A P 1A P 2B P 1B 



Rotor (shown in 

Rl ^ ^ reference position 


Reference point 1 



Rotor position, 6, deg 

(e> Step 4: salient poles Pj A and P 2A are north magnetic poles; P 1B and P 2B 
are south magnetic poles. 

Figure 5. - Resultant torque as function of rotor position lor tour-step sequence. Arrows indicate direction 
of developed torque about equilibrium points. 


In addition, the magnitude of the torque depends upon the number of phases energized, 
since the magnetic structure of the motor will be used more efficiently when all the 
phases are energized. 

For a given stepping motor, each rotor movement of 1/ 4 tooth pitch corresponds to 
an angular rotation cp in mechanical degrees. Define an angle 9 such that, when the 
rotor moves 1/4 tooth pitch, 6 varies by 90°. The angle 6 is related to the mechanical 
angle of rotation cp taken during each step by the number of teeth on the rotor. 


0 - (9) 

The resultant developed stepping motor torque from each step command can be general- 
ized by assuming that the developed torque varies sinusoidally when the rotor moves 
1/ 4 tooth pitch. By superimposing the torque produced by each pole pair, the developed 
torque can be expressed by 


T D = T max sin 6 


( 10 ) 


Figure 5 shows the resultant curves for developed torque as a function of position for 
a four-step sequence in a synchronous inductor motor. In figure 5(a), equilibrium point 
X is an arbitrary reference point. As the rotor is moved, the resultant developed 
torque is given by 


T D = " T max sin 6 


( 11 ) 


Figure 5(b) shows the resultant torque curve for the first step command. The torque is 
given by 


= - T max sin <« - 90 °> 

(12) 

t D - T max cos 6 

(13) 


Figure 5(c) shows the resultant torque curve for the second step command. The torque 
is given by 


T D = - T max sln < e - 180 °> 


T D ’ Tmax sln 6 


(14) 

(15) 


Figure 5(d) shows the resultant torque curve for the third step command. The torque is 
given by 
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T D = ~ T max sin ( 0 - 270 °) 


(16) 


T D = “ T max cos 0 


(17) 


Figure 5(e) shows the resultant torque curve for the fourth step command. The torque 
is given by 


T D = - T max sl "< e - 360 °) (18) 

T D = -T max sin 6 (19) 

This is the same result as for the reference point, equation (11). Thus, the torque de- 
veloped in a PM stepping motor for any step command sequence can be given by a com- 
bination of the four equations: 


T D " _T max sin 6 

(20) 

T D = T max cos 0 

(21) 

T D = T max sin 0 

(22) 

T D = - T max cos 0 

(23) 


Equation (7) stated that the maximum torque developed per phase was related to the 
stator current by 


T max = Kt 1 ® f 24 ' 

If the time between steps is long compared with the rise time for the stator current, the 
equation for T___ can be written 

IllaA 

T max ” K T r l = K t I 2 ■ K t I (25) 

The developed torque for one phase becomes 

T d = K t I sin e (26) 
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ment 


Figure 6. - Torque as function of rotor position for generalized permanent 
magnet stepping motor - both phases energized. 


With the aid of figure 6, equation (26) can be extended to include the case when both 
phases are energized and the magnitudes of the stator currents are assumed to be equal: 


t D = [-K t I sin(e - 45°)] + [-K t I sin(0 + 45°)] 
Tp = -K,j,l|sin(0 - 45°) + sin(0 + 45°)| 


T D = - V2 K t I sin 9 


(27) 

(28) 
(29) 


The minus sign in equation (29) is due to the reference point chosen. It indicates that 
for a clockwise (CW) 6, the developed torque will move the rotor in a counterclockwise 
(CCW) direction. 


PM Stepping Motor Single-Step Analysis 


A schematic representation of the PM stepping motor is shown in figure 7. The 
voltages supplied to the stator windings are given by 


dUt) 

E x (t) = Rl^t) + L + E y (30) 
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( 31 ) 


dl 2 (t) 

E«(t) = RI 2 (t) + L — - — + E 
1 * dt v 

The induced voltage E y generated in the stator is due to the magnetic rotor moving rela- 
tive to the stator magnetic poles. 

It is assumed that the rotor velocity vector remains perpendicular to the stator flux 
during the entire step. However, due to the geometry of the stator, the magnitude of the 
stator flux density is distributed in the stator. When this distribution is assumed to be 
sinusoidal and Maxwell’s equations are used, the induced voltage can be given by 

E = ^ sin 6(t) (32) 

v n rt dt 

Thus, the differential equation for the stator voltage per winding becomes 


E x (t) = m 1 (t) + L 
E 2 (t) = m 2 (t) + L 


+ K ae(t) gin 

(33.) 

dt v dt 


^2^ + k cos 0(t) 

(34) 

dt v dt 



In responding to a step command, the developed torque produced by the motor must 
overcome the mechanical load placed on the rotor. The mechanical load includes 

(1) The inertial torque j[d 2 <p(t)/dt^], which includes the inertia of the motor and 
the inertia of any load imposed on the stepping motor shaft 

(2) The viscous damping torque D[d<p(t)/dtj, which includes the motor damping and 
any viscous damping on the motor shaft 

(3) Torque due to Coulomb friction Tj •rjd.<^’(t)/dtj^[j dq?(t)/ dt Q | (Coulomb friction is 
defined as a friction independent of the magnitude of the relative velocity of the 
surfaces in contact, but dependent on the direction of the relative velocity. ) 

(4) External load torque applied to the motor shaft T L (t) 

In referring to figure 7, the equation of motion of the PM stepping motor can be ex- 
pressed by 


t D = K T I 2 (t) cos e(t) - Krpl^t) sin 6(t) 

de(t) 


j d^e(t) , d de(t) . ^ 

— — . ■ — 4 T 1 r 


dt 


N RT dt 2 


Nrt dt 


de(t) 


dt 


+ V-t) 


(35) 
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Figure 7. - Schematic representation of permanent magnet stepping motor for a 
single step. 


Linearized Model 

Equations (33) to (35) are nonlinear differential equations because of the torque- 
producing term in each of the equations. These equations can be linearized about an 
operating point by expanding them in a Taylor series (ref. 10) of the form 

V X 1’ X 2> = f l«l, 0> X 2, 0> + ^T (X 1’ X 2> 

Linearizing equations (33) to (36) results in 

K T I 2 (t) cos 0(t) - K T Ij(t) sin 0(t) = K T I 2 0 cos 0 o " K T I 2, 0 sin 6 0 

+ K t cos 0 q A 1 2 (t) - ILjJj 0 sin 0 Q 
- K t I 1} 0 cos e Q A0(t) - K t sin e Q Alj(t) (37) 


3fi 


UJ-i 

(X, - X. q) + — — (x. , x«) 
1 1 >' ) aXo 1 1 


L i,o 

[ 2.0 


(X 2 “ X 2, 0 ) + ' 

X l,0 

X 2. 0 


(36) 
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Transfer Function About an Operating Point at Beginning of a Step 


The linearized equations are first considered about an operating point at the begin- 
ning of a step. The constraining equation neglecting friction becomes 
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Now assume that the voltage applied to winding 1 is held constant and the motor is stepped 
by energizing winding 2. Assume the initial conditions 


(si- 


d 2 e 


— In n — ® 

2i 1,0 


<dt /o 


Also, since the voltage applied to winding 1 is being held constant, AEj(t) = 0. From the 
differential equation at the operating point 


'VljO sin 6 0 - t l,o 


(42) 


Equation (42) establishes at the beginning of the step 


0 n = sin-1 

-K T I 1, 0 


(43) 


For an excursion about the operating point, equations (37) and (38) can be combined 
as follows: 


K t cos 6 q AI 2 (t) - K t Ij 0 cos 0 O A0(t) - K t sin 0 Q Al^t) - A 

n rt 


d 2 6(t) 


Ldt 2 J 


+ ■ 


D 


N 


RT 


X A 


d0(t) 


dt 


+ AT L (t) (44) 


Making use of the identities 


sin 0 q = 


1 L, 0 
K T T 1,0 


COS 0Q = 


_V k t : m- 


L, 0 


K T I 1,0 


(45) 


(46) 


equation (44) can be written 
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N 


RT 


d 2 0(t) 


Ldt 2 J 


N. 


RT 


{%)+ (V 4 1 ?, 0-^,0) 


A0(t) = 


Kn 


V* 


2 t 2 

Th,0 


L, 0 


¥l,0 


AI«(t) 


+ ^0 AI x (t) - AT L (t) (47) 

h,o 


Also, equations (39) and (40), for an excursion about the operating point, can be written as 


AE 2 (t) = RAIg(t) + L A 


v 

RT 




o T L,o A 


¥i, 0 


0 = RAI 


[ dt J N 

«,lM M 

L dt J n rt k t x i, 0 L dt J 



(48) 


(49) 


A transfer function relating rotor position, stator voltage, and load torque can be de- 
rived by taking the Laplace transforms of equations (47) to (49). The derivations of the 
Laplace transformations are found in appendix B. The resultant linearized transfer func- 
tion about an operating point at the beginning of a step is as follows: 



(50) 


If the initial load torque T^ q is 0, equation (50) can be reduced to 

K t |aE 2 (S) + - (R + LS) AT l (S) 


A0(S) = 


JLS 3 + (DL + JR)S 2 + 


f 


/ k t K v 

[dr +-L-1 

\ n rt 


(51) 


+ LKt i i,o)S + RK t i i,o 
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Equations (50) and (51) show that the resultant linearized transfer function about an op- 
erating point at the beginning of a step results in a third-order system. The resultant 
linearized transfer function can be used to examine the effect of an initial value of load 
torque on the step response. It can also be used to examine the effects of inductance, 
generated back electromotive force (emf), inertia, and damping on the step response. 
However, the choice of an operating point at the beginning of a step makes experimental 
verification of the resulting model difficult. 


Transfer Functions About an Operating Point at End of a Step 

Consider the transfer function of the PM stepping motor about an operating point when 
the motor is in equilibrium at the end of a step with both phases energized. This opera- 
ting point is chosen because the resulting transfer function can be easily verified experi- 
mentally. At the end of the step, the initial conditions are |(dfc»)/(dt)] 0 = ~(d 2 0)/(dt 2 )J o = 0 
and 1^ q = Ig q = I. For an excursion about the operating point the constraining equation 
at the operating point is 


K T*2, 0 cos 6 0 ~ K T*1, 0 sin 0 O " T L, 0 


(52) 


or 


T 

cos e 0 - sin e 0 = -^2 

iVrpl 


(53) 


Equation (53) establishes 0 Q . If T L Q = 0, then 0 Q = 45°. The linearized equations 
about the operating point equating friction are 

-K T ljsin 0 Q + cos 0 q] A0(t) + K T cos 0 Q Alg(t) - K T sin 0 Q Alj(t) 


J A 

d 2 0(t) 

D A 

4- A 

d0(t) 

n rt 

Ldt 2 J 

n rt 

dt 


+ AT L (t) 


(54) 


AE 2 (t) = R AI 2 (t) + L A 


dl 2 (t) 

K 

+ ^ ros G A 

d0(t) 

dt 

T vUO W A 

N u 

iN RT 

dt 


(55) 
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AEj(t) = R Alj(t) + L A 


dl^t)] K y 


— sin 0 q Aj 
RT 


de(t) 


dt 


The Laplace transformations of equations (54) to (56) are 


A0(S) = 


K t cos e 0 AI 2 (S) - K t sin d Q AI 1 (S) - AT L (S) 


JS* . DS 


+ + k t I |sin 0 O + cos 0 O ] 


N RT n rt 


(56) 


(57) 


ai 2 (s) = 


E K 

AE 0 (S) + — cos 0 n S A6(S) 

S n RT 
R + LS 


(58) 


AI X (S) = 


E 1 0 K v 

AE^S) + -^-2-- — 5L sin 0 n S A6(S) 

s n rt 

R + LS 


(59) 


Combining equations (57) to (59) gives the resultant linearized equation for both phases 
energized about an operating point at the end of a step 


K t cos 


A0(S) = 


R 


cos e n r E„ n “| K™ sin 6L E., n 

2 AE,(S) +— til - 2 AE-(S) +— L2 - AT t (S) 

+ ls L sj r + ls L 1 s J g 

DS | K T K v 
RT N RT^ R + LS ^ 


C + iz H 2 e ° - s1 " 2 e °> + v( sln + cos e °) 


(60) 


For no initial load torque (T L Q = 0), the constraining equation (eq. (53)) yields 

A * 


0 Q = 45 . Thus, 


0. 707K n 


A0(S) = 


R + LS 


E, 


AE 9 (S)+^M 


0. 707K n 


R + LS 


AE 1 (S) +_!i2 

1 c 


- AT l (S) 


(61) 


JS1 + _DS_ + V^k t I 

n rt n rt 
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For load torque disturbances about the operating point at the end of the step and the initial 
load torque equal to 0, the transfer function is a second-order system 

-- at l (s)(n rt ) 

0(S) = 2 — — (62) 

2 D , V^T^RT 

JS J 

This result is important for two reasons. First, it allows the stepping motor to be de- 
scribed in terms of a natural frequency and damping ratio. Second, these parameters 
can easily be determined experimentally. 

Equation (62) has been developed in terms of 0, such that one step equals 90°. 
Equation (62) can be written in terms of the actual mechanical degrees for a particular 
stepping motor. Since 


d 2 <?(t) [ D d<p(t) , T 1 
dt 2 J dt J 


e(t) = 


AT L (t) 

J 


(63) 


it follows that 


- AT l (S) 

A^(S) = (64) 

c 2 D , N RT ^ K t J 
JS J 

The natural frequency obtained from equation (64) is given by 



(65) 


The natural frequency of the PM stepping motor is increased by 

(1) Decreasing the inertia 

(2) Increasing the number of rotor teeth 

(3) Increasing the magnitude of the stator current 

The upper bound on the natural frequency due to inertia is often limited by the load iner- 
tia. However, minimizing the rotor inertia increases the natural frequency when the load 


20 


inertia is not the limiting factor. The natural frequency of the stepping motor can also be 
expressed in terms of the stall torque Tg with both phases energized 

T g = V2K t I (66) 


Thus, 

\ 


.. - RT T S 

“*-y j 

(67) 

The damping ratio of equation (64) is given by 


N = - 
w J 

(68) 

^ _ D 

(69) 


2 V jn rt t s 


The damping ratio is increased by increasing the viscous damping, or by decreasing the 
inertia. The damping ratio is decreased by increasing the number of rotor teeth or by 
increasing the maximum torque by increasing the stator current. 

The natural frequency and the damping ratio of the linearized model of the PM step- 
ping motor can also be evaluated by considering a small change in A6 from the equilib- 
rium point. Consider the linearized equations about the point of stable equilibrium with 
the stator current held constant and no applied load torque. The linearized equations re- 
duce to 


-K T I(sin 9q + cos 0 Q ) A0(t) 


J A 

d 2 0(t) 

+ D A 

de(t) 

n rt 

_dt 2 _ 

n rt 

dt 


(70) 


The Laplace transform of equation (70) is 
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-K T I(sin 0 q + cos 0 q) A0(S) = ^ 

n rt 



+ 


D 


N 


RT 


S A0(S) - A0(O) - 



( 71 ) 


For a small initial displacement A0(O) and initial conditions 6q = 45° and 
A[d0(O)/dt] = (d0/dt) Q = (d 2 0/dt 2 ) o = 0, equation (71) can be rewritten as 

_N RT ^2 K T X A6 ^ = ^ Js2 + DS ^ ” ( JS + A0 (°) 


(72) 


where 


A0(S) = 


S + 


?) 


A0(O) 


c 2 D , V2 K T m R T 
JS J 


(73) 


Equation (73) can be used to evaluate the natural frequency and the damping ratio of the 
linearized model about the stable equilibrium point. 


PHASE PLANE ANALYSIS 

The linearized mathematical model of the PM stepping motor was developed to ana- 
lyze the single-step response of the PM stepping motor. For multistep operation, the 
response of the PM stepping motor for a given load torque is governed by 

(1) The rate of the input pulses 

(2) The current transients in the stator windings 

(3) The mechanical parameters of the rotor 

The linearized model cannot be used to analyze multistep responses because the model 
fails to account for motor failures caused by excessive input stepping command. 
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The transient response of the stator current has a significant effect on the stepping 
motor response as the stepping rate is increased. Even if the back emf induced in the 
stator is not significant, the current turnon transient affects the maximum torque de- 
veloped by the motor. Consider the electrical (L/R) time constant of the stator windings. 
If the time between the application of step commands approaches this time, the current 
will not reach its expected value. The maximum torque developed by the motor is re- 
duced, and thus, the natural frequency of the stepping motor is reduced. The amount of 
load torque which the motor can step against is also reduced. 

In general, stepping motor drive circuits are designed to compensate for stator cur- 
rent transients. The drive circuit tends to act as a constant current source by control- 
ling the stator current independently of the inductance or generated emf. A typical con- 
stant current drivS source used to compensate a stepping motor is presented by Zeller 
(ref. 11). 

For multistep analysis, it is assumed that a drive circuit is used which approximates 
a current source for stepping motors. The stator currents are assumed to be constant 
and to have magnitudes equal to I or -I during each step. 


Analytical Model Assuming Constant Current Source 

The assumption of a constant current source simplifies the analytical model of the 
PM stepping motor. It also permits the response of the motor to be studied for a series 
of input pulses. Because the current is independent of dl(t)/dt and d0(t)/dt, the stator 
winding equations become 


Ej = ±IR 

(74) 

E 2 = ±IR 

(75) 


The sign in each equation is determined by the input pulse sequence. 

The differential equation of rotor motion can be reduced to a second-order nonlinear 
differential equation. For the rotor at the equilibrium position, the developed motor 
torque has been shown to be 

T d = K t I sin[e(t) - 45°] + K T I sin[fl(t) + 45°] (76) 


A step command requires that the rotor angle 9 advance 90°. The developed torque 
generated by the motor to accomplish this step is given by 
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or 


T d = K t I sin 1 90° - [fl(t) - 45°]) + K T I sin ^90° - [(9(t) + 45°]} 


(77) 



Equation (78) reduces to 

T d = V2K t I cos 0(t) (79) 

Equating equation (79) to the mechanical load of the motor results in 


V2K t I cos 6( t) 


J d 2 g(t) D 
N RT dt 2 N RT 


d0(t) , 


dt 


dg(t) 

dt 

de(t) 
dt ; 


+ t l W 


(80) 



Figure 8. - Schematic representation of rotor command positions. Multistep operation. 
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Multistep Analysis Assuming Constant Current Source 

Equation (80) can be extended to include multistep operation. Each step command 
calls for the rotor equilibrium position to shift by 90°. Examine figure 8 and assume the 
rotor is initially at 0°. The equations for the developed torque for a series of steps are 


then 

Step 1: T d = K t I sin j 90° - [e(t) - 450} + K T I sin j 90° - (j9(t) + 45°]} (81) 

Step 2: T d = K T I sin |l80° - [fl(t) - 450} + K T I sin j 180° - [e(t) + 450} (82) 

Step 3: T d = K T I sin |270° - [fl(t) - 45*0} + K T I sin | 270° - [e(t) + 450 } (83) 

Step 4: T d = K T I sin |360° - [e(t) - 450} + K T I sin j 360° - [e(t) + 450} (84) 

Step 5: T d = K T I sin ^450° - [e(t) - 450} + K T I sin ^450° - [e(t) + 450} (85) 

and so forth. These equations reduce to 

Step 1: T d = V2 K t I cos 0(t) (86) 

Step 2- T d = V2 K t I sin 0(t) (87) 

Step 3: T D = - V2 K T I cos 6(t) (88) 

Step 4: t d = - V 2 K T I sin 0(t) (89) 

Step 5: T D = V2 K T I cos 6(t) (90) 

and so forth. The basic equation for multistep operation becomes 
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V2K t I sin 


[«, - «»>] J 


d 2 fl(t) 


N 


RT dt 


N 


RT 


dfl(t) 

dt 


+ T x 


d6>(t) 

dt 

dO(t) 

dt 


+ T t (t) 


(91) 


Normalized Model for Multistep Operation 

The second-order nonlinear differential equation obtained in equation (91) lends itself 
to phase plane analysis. This equation can be considered as the transient solution to a 
set of initial conditions. By specifying the initial conditions 0(0) and d0(O)/dt, the solu- 
tion for all positive time is completely determined. An advantage of the phase plane 
analysis is that trajectories representing the transient solution of the equation can be de- 
termined without solving the equation for the dependence of position on time. 

Equation (91) can be normalized to reduce the number of parameters. By using the 
linearized analysis, the normalization can be performed to allow the use of experimental 
data. Dividing equation (91) by the stall torque ^T g = yfl K^l) results in 



J d 2 0(t) | D 

N rt V 2K t I dt 2 N rt V 2K t I 


d0(t) , T f 

d0(t) 

dt 

t l 

1 

dt V 2 k t i 

d0(t) 

dt 

T 

V2K t I 


(92) 


Recall that the natural frequency of the linearized second-order model is given by 

1 /N r V2K I 

p 

Equation (92) can be simplified by defining a new time parameter in terms of the natural 
frequency: 


T = 


V 


(94) 


d _ 1 d _ / J d 

d T u> N dt y yfe k t I dt 


(95) 
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dr 2 N rt V 2K t I dt 2 


( 96 ) 


D 


D 


n rt V^Rt 1 dt ^/n rt V2K t IJ 


A = D — 
dr dT 


(97) 


Equation (97) can also be defined in terms of the damping ratio 


?=» 


D 


„ V N R T _ D 


2o^.J 




RT V^K t I 


(98) 


or 


D = 2? 


(99) 


Also, normalizing the friction and load torques gives 


T f = 


T f T f 


-\[2 r:r T T 


and 


T L = 


K t I x s 


t l t l 


V 2 K t I t s 


( 100 ) 


( 101 ) 


Thus, equation (92) can be written with 6 and 8 c specified in radians 

d0 

sin(6 - 0) = ^ + D — + T- -^—+ T t 
c dj 2 dr f I d0 1 L 


dT 


( 102 ) 


Time can be eliminated from equation (102) by defining 
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v = de 

dr 


(103) 


> 


dV = dV d0 = y dV 
dT d0 dT d 0 


Rewriting equation (103) in terms of position 0 and velocity V results in 

sin(0 - 0) = V — + DV + T f + T t (104) 

c d0 1 |V| L 

The phase plane solution to equation (104) is a plot of velocity V as a function of po- 
sition 0. The initial conditions V(0) and 0(0) locate an initial point in the phase plane. 
The trajectory through this point describes the response of equation (104) for all positive 
time. The slope of the trajectory is given by 


dV 

d0 


sin(0 c 


- 0) - DV - T 


V 



(105) 


For a given motor with a specified value of friction, the slope of the trajectories is a 
function of three variables, 


^ = f(0,V,T T ) (106) 

d0 ' Lj/ 

Thus, the unique response to a set of initial conditions must be specified in terms of the 
load torque, assuming that the motor is loaded by a speed independent torque of magni- 
tude T l . In the normalized equations, 0 and 0 c must be specified in terms of radians. 


Phase Plane Analysis by Delta Method 

There are several graphical methods to determine the phase plane trajectories. For 
single-step analysis, the most convenient method is the delta method (ref. 10). The delta 
method assumes that the trajectories can be approximated by arcs of circles with origins 
on the 0-axis. In figure 9(a), 5 gives the center of a circular arc through point Pj. 
Comparing similar triangles, 


28 


dV _ 0 _ 0 + 5 


(107) 


de v v 

-6 = 0 + V — (108) 

de 

To determine a trajectory using equation (108), first select a point P^ and solve for 
delta. The value of delta will determine the center of an arc r c that will pass through 
point Pj. Notice that the center of the arc always lies on the 0-axis. After drawing a 
segment through point Pj, select a point P2 on the segment. Compute a new delta and 
new arc center, etc. Figure 9(b) illustrates the method. Using equation (104) with 


v 



(a) Determination of center of circular arc through point Pj. 



0 45 90 

Rotor position, deg 

(b) Use of delta method. 

V dV/de - cos(9 - V). 

Figure 9. - Graphical solution of a phase plane trajectory by delta method. 
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V(0) = 0(0) = Tj = = 0, 6 c = v/2 , and D = 1. 0 results in 


= cos 0 - V (109) 

d0 


-6 = 0 + V — — = 0 + cos 0 - V (110) 

d 0 

The process for computing the trajectory continues until the equilibrium point is reached. 

Time can be computed along the trajectory (ref. 12). The relation between time, po- 
sition, and velocity is given by 


T = 



d0 


(HI) 


For a segment along the trajectory, 


t N+1 ' 



( 112 ) 


or 


t N+1 


%+l ~ ^jST 

|( v n + i +v n) 


(113) 


The total time is the summation of the increments of time along the entire trajectory. 

Single-step response. - For a single-step response the phase plane is useful in ana- 
lyzing 

(1) The transient response to a set of initial conditions 

(2) The effect of load torque T L on the transient response 

(3) Other nonlinear effects, in particular Coulomb friction, provided they are built 

into the model 

The equation for a single-step response was given in equation (104) as 

sin(0 - 0) = V — + DV + T f — + T t (114) 

c d0 1 | V | L 

The step response will first be analyzed for T^ = T^ = 0. Equation (114) reduces to 
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(115) 


sin(e - e) = V— + DY 
c de 


Assume that the rotor is initially at rest and that the step command will move the rotor 
from 0° to 90° as shown in figure 8. The equation for the step becomes 




+ DV 


(116) 


or 


cos 0 = V — + DV (117) 

de 

For a given value of the normalized damping ratio D, equation (117) describes a trajec- 
tory in the phase plane for a given set of initial conditions. 



0 45 90 135 

Rotor position, 9, deg 


Figure 10. - Effect of normalized damping on single-step response. Nor- 
malized frictional torque Tj, 0; rotor command position 0 C , 1.57 radians. 
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Figure 10 shows the effect of D on the step response. Increasing D makes the re- 
sponse less oscillatory and increases the rise time for each step response. In figure 10, 
the normalized time is given for each step response when 6 reaches 1. 50 radians or 
about 95 percent of the command value of 1. 57 radians (90°). The normalized time var- 
ies from 1. 97 to 5. 74 as D varies from 0. 25 to 2. 0. 

Effect of load torque . - The addition of load torque causes an offset in the equilib- 
rium position of the rotor. Load torque also slows down the step response. With load 
torque applied, the initial position of the rotor is determined by the equilibrium point of 
the previous step. Figure 11 illustrates this effect. Figure 11(a) shows the developed 
motor torque when T L = 0. Assume the equilibrium position with no load torque before 
the step command is given is 0°. The previous step command required to reach this po- 
sition is given by 


sin(0° - 6) = V — + DV (118) 

d6 

-sin 6 - V + DV (119) 

de 


The equilibrium position is given by 

- sin 6=0 j 


or 


> 


6 = 0°(0 radians) 


(120) 


These form the initial conditions for the next step. Notice that at this position no motor 
torque is developed. If the rotor is moved about this point, the motor develops a torque 
which tends to drive it back to the point of equilibrium. When the next step is com- 
manded, maximum motor torque is applied to the rotor for acceleration. As the rotor 
moves toward 90°, the amount of available torque is reduced until the 90° position is 
reached. At this point, the developed torque is again 0. 

Figure 11(b) shows the developed motor torque for a normalized load torque of 0. 4. 
For the previous step, the equilibrium position has shifted 23. 6°. This shift is due to the 
fact that some motor torque is required to hold the load torque. The initial equilibrium 
position of -23.6° was determined from 

-sin 6 = V — + DV + 0. 4 (121) 

d6 
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(a) T l = 0. 



(b) T L = 0. 4. 

No torque available 



(c) Tl = 0.707. 

Figure 11. - Effect of normalized load torque T L on developed 
motor torque for a step command. Arrows indicate direction 
in which developed motor torque acts. 



When the rotor reaches equilibrium, 


-sin 6 = 0. 4 

6 = -23.6°(0.645 radians) 


> 

j 


(122) 


The shifted equilibrium position reduces the amount of motor torque that is available to 
accelerate the rotor to the next step. The torque that is available for acceleration is 
given by 


cos 0 - T^ 



+ DV 


(123) 


The amount of motor torque available for acceleration decreases with increasing load 
torque until T^ = 0. 707. From figure 11(c) for a load torque of 0. 707, the rotor position 
offset is 45°. When the next step is commanded, 


cosf- - 0. 707) = V — + DV (124) 

\4 / d6 


0 = V — + DV (125) 

de 

there is no motor torque available to accelerate the rotor, and a step command cannot 
be accomplished. The motor can statically hold a load torque equal in value to the stall 
torque of the motor. But the motor cannot step a load torque that is equal to or greater 
than 0. 707 of the motor's stall torque Tg. Thus, a normalized load torque of 0. 707 will 
cause the PM stepping motor to fail to respond to a step command. 

Figure 12 shows the step response in the phase plane for various values of load 
torque. Normalized damping I) is given as 1. 0. Notice that although the initial and final 
rotor positions are offset by the load torque, the distance between these positions is 90° 
(1. 57 radians). Thus, load torque does not alter the size of the step, just the initial and 
final positions of the rotor. 

Figure 12 illustrates the effect of load torque on the speed of response. The maxi- 
mum velocity reached is reduced with increasing load torque. This is a direct result of 
the available motor torque to accelerate the rotor. Notice that with increasing load 
torque the amount of position overshoot is reduced. Also shown in figure 12 is the nor- 
malized time for the rotor to move 50 percent of its step. The time is increased with in- 
creasing load torque. 
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Normalized rotor velocity, V 


B (1. 33) -a 

Normalized 

time, 

T 


\ 

\ 

< AJjL5B^ 

Normalized 
load torque, 

V T. 


Rotor positron, rad 


0 45 

Rotor position, deg 


Figure 12. - Effect of normalized load torque on single-step response. Normalized damp 
ing, 1.0; frictional torque, 0. 


Normalized 

frictional 


torque, 


Rotor positron, 8, rad 


45 90 135 

Rotor position, 8, deg 

Figure 13. - Effect of normalized Coulomb friction on single-step response 
without load torque. Normalized damping, 0. 5; normalized load torque, 0. 





Effect of Coulomb friction. - Coulomb friction changes agiven step size. Figure 13 
illustrates the step response for various values of T f with Tj = 0 and^ D = 0. 5. Near 
the command position, the developed motor torque becomes small. If Tj > cos 8 when 
V = 0, the friction will hold the rotor, causing an offset from the command position. 
Thus, the equilibrium position lies in a band about the command position. The size of 
the band is determined by the amount of friction present. The speed of response is also 
reduced by friction. The reduction is small, however, unless a large amount of friction 
is present. 



0 45 90 

Rotor position, 0, deg 


Figure 14. - Effect of normalized Coulomb friction on single-step response 
with load torque. Normalized damping, 0.2; normalized load torque, 0.5. 


Combined effects . - Figure 14 illustrates the effect of friction combined with load 
torque. The effects are similar to those experienced with load torque alone. The final 
rotor position for a given step is determined by the values of load torque and friction 
present. The effect of Tj on stepping motor operation will not be considered further in 
this report since the frictional torque magnitude is generally small and its effects are 
somewhat similar to the load torque effects. 

Multistep analysis . - The real advantage of using the phase plane in analyzing the 
stepping motor is that the analysis can be extended to include multistep operation. In 
multistep operation, each input step to the motor commands the rotor to move to a new 
position 90° away from the existing command position. The dynamic equation that re- 
lates the motor movement is equation (104). If the time between the application of input 
steps is long, the rotor may reach the new point of equilibrium before the next step com- 
mand is applied. As the rate of step commands is increased, the rotor may not reach 
the new equilibrium point. At the instant the next step command is applied, the position 
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and velocity of the rotor form initial conditions for the new equation of rotor motion. The 
step commands are applicable to unidirectional or bidirectional rotor movements . 

The phase plane for multistep analysis can be considered as a series of trajectories 
for each possible dynamic equation given by the step commands. Time is computed along 
the trajectories and is used to indicate the time for the application of the next step com- 
mand. As an illustration, consider the following series of input commands with Tj = 0: 


Step 1: 

sln (H 

l = cos 8 = V — + DV + T t 
1 d0 L 

(126) 

Step 2: 

sin(7r - 8) = sin 8 = V — + DV + T t 

d6 ^ 

(127) 

Step 3: 

sin(~ - 8 s 
\2 J 

1 = -cos 8 = V — + DV + T t 
' de L 

(128) 

Step 4: 

sin(27r - 8) = -sin 6 = V — + DV + T t 

d6 h 

(129) 

Step 5: 

sin(— - £ 

V 2 

n = cos 8 = V — + DV + T t 
/ d8 ^ 

(130) 


Figure 15 shows a phase plane portrait for this series of steps. The trajectory is for 
D = 0.25 and T^ = 0. Consider a normalized stepping period of At = 1. 31. 

The step commands are applied at the following normalized times: 


Step 1: 

r = 0, 6 c = 90°^- radiansj 

Step 2: 

r = 1. 31, 8 c = 180° (v radians) 

Step 3: 

t = 2.62, & c = 270°^— radians^ 

Step 4: 

t = 3. 93, 8 c = 360° (2n radians) 

Step 5: 

t = 5. 24, & c = 450°^— radians^ 
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Normalized rotor velocity, V 
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Command position: First step (90°) Second step U80°) Third step (270°) Fourth step (360°) Fifth 

l i i i i i i ( 4 f° 

1 2 3 4 5 6 7 8 

Rotor position, 0, rad 

Figure 15. - Phase plane portrait for multistep input with no load torque. 
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Command First step 
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Second step 
(180°) 


Third step 
(270°) 


Fourth step 
(360°) 


0 1 2 3 4 5 6 7 

Rotor position, 8, rad 

Figure 16. - Phase plane portrait for multistep input with normalized load torque equal to 0. 2. 



The initial conditions for step 1 are V(0) = 0(0) = 0. The trajectory is governed by 
equation (126) and is shown in figure 15 between 0=0, V = 0, and point A. At point A, 
the total response time r is 1. 31, and step 2 is applied. The values of V and 0 at 
point A become the initial conditions for the response of step 2. The equation governing 
step 2 is equation (127). The trajectory moves from point A to point B. At point B, the 
total time r is 2.62. The third step command is applied and is governed by equa- 
tion (128). The values of V and 0 at point B become the initial conditions for equa- 
tion (129). This process is continued for each step applied. The rate of stepping can be 
changed by simply changing the time between steps At. 

The dashed trajectories in figure 15 show the phase plane trajectories if the suc- 
ceeding commands are not applied. Curve 0A' is for a single-step command, OB' is for 
a two-step command, and so forth. The response is oscillatory in each case, but each 
trajectory reaches its correct command position. The oscillations of the trajectories 
as they approach the equilibrium position are a function of the maximum velocity reached 
during the stepping sequence. 

Figure 16 is a phase plane portrait of the same stepping sequence but with a load 
torque T^ of 0. 2. The addition of load torque causes the rotor velocity to build up at a 
slower rate than in the case of no load torque. Without load torque, the rotor velocity 
peaks during the third step. With the load torque, the velocity does not peak until the 
fifth step. The dashed trajectories in figure 16 show the phase plane trajectories if the 
succeeding commands are not applied. A comparison of these trajectories clearly shows 
the applied load torque slows down the speed of response. 


Phase Plane Analysis by Method of Isoclines 


The phase portrait for the differential equation of the PM stepping motor can be used 
to determine whether the stepping motor will fail for a given stepping rate with a given 
value of load torque. The delta method provides a convenient way of constructing the 
phase trajectories for a given set of initial conditions and lends itself to computer solu- 
tion. However, the entire family of phase trajectories can be graphically examined by 
developing the phase portrait using the method of isoclines (ref. 12). This method de- 
velops the phase plane by plotting curves, called isoclines, of constant trajectory slope. 

The slope of the stepping motor differential equation is given by 


dv _ sin(0 c - 0)- T l - 
d0 V 


(131) 


The slope of the differential equation must be specified in terms of a given D and T^. 
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The value of D is determined by the specific system in which it is applied. For the fol- 
lowing analysis, D is chosen as 0. 25. The load torque will be neglected in the following 
analysis. For a step command of 90°, the slope is given by 

— = £2±JL- 0.25 (132) 

de v 

Figure 17 shows a plot of the isoclines for equation (132). A phase trajectory, for a 
given set of initial conditions, is plotted by connecting short line segments with the proper 
slope at intervals along the isoclines. A trajectory with initial conditions 6(0) = 0 and 
V(0) = 0 is plotted in figure 17. 

Figure 17 shows the existence of singular points in the phase plane. A singular point 
is a representative point (9q, Vq) in the phase plane for which (d0/dt) = (dV/dt) = 0 
(ref. 12). For a singular point, the direction of its tangent is indeterminate, and its tra- 
jectory degenerates into the singular point itself. From figure 17, the singular points 



Rotor position, 0, rad 

I 1 I I _J 

-90 0 90 180 270 

Rotor position, 0, deg 

Figure 17. - Isoclines for dWd0 = cos 0 - 0.25. 


V 


V 




(a) Focal point. Indicates (b) Saddle point. Indicates 

singular point of stable singular point of un- 



-270 -180 -90 0 90 180 270 360 450 540 630 

Rotor position, 0, deg 


(c) Singular point locations for stepping motor equation. 

Figure 18. - Types and locations of singular points for normalized stepping motor 
equation. 


are found, for example, at 0 = ±90°, ±270°, and ±450°. The singular points for these posi- 
tions are either focal points or saddle points (ref. 13). Focal points are singularities in 
which trajectories approach in a decreasing spiral manner, as shown in figure 18(a). 

Focal points form points of stable equilibrium for the stepping motor. Saddle points are 
singularities which are approached by trajectories forming distinct sectioned curves, as 
shown in figure 18(b). Saddle points form points of unstable equilibrium. Figure 18(c) 
shows the alternating locations of focal points and saddle points. The focal points are 
360° apart and 180° from the saddle points. 

Multistep analysis . - Multistep operation can be analyzed for a given stepping se- 
quence by plotting the isoclines of the resulting equations for the desired step commands. 
For a given stepping sequence, set of initial conditions, and stepping rate, the phase tra- 
jectory is developed by selecting the set of isoclines corresponding to the given step com- 
mand. The phase trajectory is then drawn from the initial conditions. Time is computed 
incrementally along the trajectory. When the computed time equals the stepping period 
for the given command step, the next step command in the sequence is applied. The point 
on the trajectory at which the computed time equals the stepping period becomes the ini- 
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tial condition for the next step command. The trajectory for this step command is plotted 
by using its corresponding set of isoclines. Time is computed along the trajectory and 
compared with the stepping period to determine the application of the succeeding step 
command. This method is continued for the given stepping sequence. 

Effect of steppi ng rate. - Figures 19(a) to (d) show the effect of stepping rate on the 
ability of the PM stepping motor to follow a set of sequential step commands. These fig- 
ures show phase plane plots developed by the method of isoclines and consider the set of 
sequential step commands that were developed in equations (126) to (130). The equations 
for the slope of the trajectories with D = 0.25 and T^ = 0 are 


Step 1: 6 = 90°(— radians\; — = - cos - - 0.25 (133) 

c \2 / d6 V 


Step 2: 6 = 180 o (t radians); — = - s * n - - 0.25 (134) 

c d6 V 


Step 3: 6 = 270°/— radians'); — = - £ os 1 - 0.25 (135) 

c \2 ) d8 V 


Step 4: 8 = 360°(2tt radians); ^ = _sin_e _ 25 (136) 

c de V 


Step 5: Q = 450°/— radians') ; — ^ - os 6 - 0. 25 (137) 

c \2 / d6 V 


For the sequential step commands, again consider a normalized stepping period of 
At = 1. 31. Recall that the step commands are applied at the following normalized times: 


Step 1: 

r = 0; 6 q = 90°^- radiansj 

Step 2: 

t = 1. 31; 6> c = 180°(7 t radians) 

Step 3: 

r = 2. 62; e Q = 270°^ radiansj 

Step 4: 

t = 3. 73; 6 q = 360°(27t radians) 

Step 5: 

r = 5. 24; 6 c = 450°^— radiansj 


44 



Step 1 calls for a position command of 90°. The trajectory for the response to this 
command is shown in figure 19(a). The rotor is assumed to be initially in equilibrium at 
6(0) = 0 and V(0) = 0. The trajectory is drawn from 6(0) and V(0) to point A. At 
point A, the normalized time r is 1. 31, and the command for step 2 is applied. Step 2 
calls for a position command of 180°. The trajectory for the response to this command 
is shown in figure 19(b) from point A to point B. At point B the computed time is 2.62, 
and the command for step 3 is applied. Step 3 calls for a position command of 270°. 

The trajectory for the response to this command is shown in figure 19(c) from point B to 
point C. At point C the computed time is 3. 93, and the command for step 4 is applied. 
Step 4 calls for a position command of 360°. The trajectory for the response to this com- 
mand is shown in figure 19(d) from point C to point D. At point D the computed time is 
5. 24, and the command for step 5 is applied. Step 5 calls for a position command of 
450°. The slope of the trajectory in response to this command is given by equation (137). 
This equation gives the same result as for the slope of the trajectory in response to the 
90° step command given by equation (133). Thus, the response for the 450° position com- 
mand can be shown in figure 19(a) from point D to point E. This procedure can be con- 
tinued for plotting the trajectory in response to any additional step commands desired. 

Failure t o execute commands . - Increasing the stepping rate decreases the amount 
of time in which the motor can respond to the step command. As a result, the rotor 
tends to lag behind the command position until sufficient velocity is built up to allow the 
rotor to reach the command position. The rotor can dynamically lag behind the command 
position by as much as 180° (two steps) and still follow the command. But once the rotor 
falls behind by more than 180°, the motor cannot reach the correct command position. 

As an illustration, consider the previous step sequence with the normalized stepping 
period At reduced from 1. 31 to 0.92. The step commands are applied at the following 


times: 



Step l': 

r= 0; 

& c = 90°^— radians^ 

Step 2’: 

t= 0.92; 

6 c = 180 °(tt radians) 

Step 3’: 

t = 1. 84; 

6 C = 270° radians^ 

Step 4’: 

r = 2. 76; 

& c = 360°(2jt radians) 

Step 5* : 

r = 3. 68; 

6 q = 450°^— radians^ 
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The response to step 1* is shown in figure 20(a) from V(0) = 0 and 0(0) to point A’. 

At point A T , r = 0. 92, and the 180° step command is applied. Notice that the value of 
velocity and position reached is considerably lower than in the previous case for t = 1 . 31. 
Figure 20(b) shows the trajectory for the response to the 180° step command between 
points A’ and B*. Note in figure 20(b) that the rotor position at point B' is lagging the 
command position by more than 90°. At point B’, t = 1. 84, and the step command for 
270° is applied. Figure 20(c) shows the trajectory for the response to the 270° step com- 
mand, between points B’ and C'. The velocity of the rotor actually decreases during this 
step, and the rotor lags further behind the command position. At point C', r = 2. 76, and 
the 360° step command is applied. Note by examining the isoclines, that had the 360° 
step command not been applied, the rotor could still have reached the correct command 
position of 270°. Figure 20(d) shows the trajectory for the response to the 360° step com- 
mand between points C' and D\ Again during the step the rotor velocity decreases. At 
point D' the rotor is approximately 180° behind the command position. However, the 
rotor can still reach the proper command position of 360°. At point D’, r = 3.68, and the 
450° step command is applied. Figure 20(a) shows the response to the 450° step com- 
mand between points D' and E’. The rotor is more than 180° behind the command posi- 
tion during the entire step. Clearly, the motor cannot reach the desired command posi- 
tion of 450°. If no further steps are applied, the rotor will come to rest at 0 = 90°, or 
four steps behind the desired command position. 

In general, the rotor will fail to execute the step command if the stepping motor tra- 
jectory crosses the position axis before the desired command position is reached. This 
can only occur if the position of the rotor lags the step command by more than 180°. The 
failure analysis must be carried out by using enough steps to show that either the motor 
fails to follow the desired command or that the motor reaches a steady- state condition 
for the given command rate. 

If a motor will follow step commands at a given stepping rate, it will follow the step 
commands regardless of the number of steps applied. Also, if it cannot follow a step 
command, the motor will continue to lose steps. The steps that it loses can be applied 
after the step commands have been completed. This can be accomplished by feedback but 
introduces more complications (ref. 1). 

The method of isoclines is not limited to analyzing step commands in one direction or 
to a constant time between steps. Any step sequence or stepping rate can be considered 
by selecting the proper set of isoclines corresponding to the step sequence. In general 
any step sequence or stepping rate can be considered by selecting the proper set of iso- 
clines corresponding to the step sequence. Thus, any input step sequence can be inves- 
tigated by this method. 
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DIMENSIONLESS CURVES FOR EXPRESSING PM STEPPING 
MOTOR PERFORMANCE 


Multistep operation in the phase plane has been shown to be a function of three di- 
mensionless parameters: 

(1) Motor damping parameter, D = 2? (138) 


(2) Time parameter 




RT ^ K T* 




(3) Load torque parameter, T^ = 


T L t l 


V 2 K t I t s 


(139) 


(140) 


Using these parameters, multistep analysis in the phase plane can be used to specify the 
maximum stepping rate for any given PM stepping motor and desired value of load torque. 

The maximum stepping rate is defined as the rate in which the motor can step with a 
given value of load torque and not miss a single step. The PM stepping motor is assumed 
to be initially at rest at a point of stable equilibrium for the given value of load torque. 

By determining the maximum stepping rate for various values of D and T^, a set of di- 
mensionless curves can be developed to predict PM stepping motor performance. 

A computer program was developed which uses the delta method to map the maxi- 
mum stepping rate as a function of D and and is described in appendix C. In the 
computer program, the maximum stepping rate is determined by an iteration process. 

A normalized stepping rate 1/Ar is chosen and expressed in terms of the time between 
application of steps At. The motor is then stepped at the given rate for a number of 
steps. The phase portrait is examined to see if the motor follows the command without 
missing a step. If the motor follows the command, a higher stepping rate is chosen, and 
the procedure is repeated until a stepping rate is found in which the motor cannot follow. 

Figure 21 shows the resulting dimensionless curves of PM stepping motor perfor- 
mance derived from the computer program. The curves show that, for a given value of 
load torque, the maximum stepping rate of a given PM stepping motor decreases as the 
damping ratio increases. This has a direct bearing on the application of the PM stepping 
motor to system design. The motion of the PM stepping motor for a single step is highly 
underdamped, yielding low values of damping ratio. The motor can be compensated to 
obtain a higher damping ratio and a better single-step response. But raising the value of 
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0 1 2 3_ 4 

Normalized motor damping, D 

Figure 21. - Dimensionless parameter curves for 
permanent magnet stepping motor. 


the damping ratio decreases the maximum stepping rate that can be obtained. Thus, for 
multistep operation, compensating for the single-step response may not be desirable. 

The phase plane analysis demonstrated that the ability of a motor to follow a set of 
step commands was, in part, determined by the position of the rotor corresponding to the 
command position during each succeeding step. If the rotor position falls further behind 
with each succeeding step, the motor will eventually miss a step. Thus, the rise time 
and, hence, the damping ratio of the rotor during each step affects the maximum stepping 
rate of the PM stepping motor. The damping ratio is given by 

C = (141) 

2Jw n 

Thus, the damping ratio can be lowered by decreasing the viscous damping or increasing 
the motor inertia or the natural frequency. However, increasing the inertia decreases 
the natural frequency since 
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(142) 


w n 


N 


r 

t 



V 2 k t i 

J 


If the natural frequency of the motor is increased, the time between steps can be de- 
creased, and the maximum stepping rate can be increased since 


Ar = Ci^At 


(143) 


At At 


(144) 


Thus, for a given load torque, maximum stepping rate is obtained by minimizing the 
damping ratio and maximizing the natural frequency. 

Figure 21 also shows that, for a given damping ratio, increasing the load torque de- 
creases the maximum stepping rate that can be achieved. Since the dimensionless curves 
assume a constant current source, the maximum stepping rate will be decreased if the 
motor is driven from a voltage source. The stepping rate in this case is decreased as a 
result of the stator winding L/R time constant and induced emf, both of which cause the 
rotor position to further lag behind the intended command position. 



Figure 22. - Experimental setup for verification of analytical model. 


51 












Normalized maximum stepping rate, 1/At 



Normalized motor damping, D ■ 2£ 


Figure 23. - Experimental verification of dimensionless 
parameter curves for normalized applied load torque 
T l of 0. 



Figure 24. - Experimental verification of normalized step- 
ping rate for various values of load torque with normal- 
ized damping D of 0.9. 


EXPERIMENTAL VERIFICATION OF DIMENSIONLESS CURVES 


Some experimental measurements were taken to check the validity of the dimension- 
less curves developed for the PM stepping motor. Figure 22 is a schematic diagram of 
the apparatus used in the experimental work. The measurements were taken on syn- 
chronous inductor stepping motors. The stepping motor drive circuit used approximated 
a current source. Damping torque was obtained by adding weights to the motor’s output 
shaft. Additional damping torque was obtained by coupling to the stepping motor output 
shaft a dc servomotor whose field windings were open- circuited. Load torque was ob- 
tained by using the servomotor in the conventional manner. 

Figure 23 shows a comparison of the experimental values obtained to the analytically 
derived curve for T^ = 0. The experimental data show reasonable agreement with the 
analytical curve. 

Figure 24 shows the comparison of experimental data with the corresponding ana- 
lytical data from the dimensionless curves for values of T^ of 0, 0.2, 0.4, and 0.6. 

The normalized damping was held at 0.9. The results show good agreement. 

CONCLUSIONS 

An analysis of the permanent magnet (PM) stepping motor can be used to study the 
effects of inertia, damping, current transients, etc. , on the stepping motor response to 
a single-step command. A nonlinear analysis which uses phase plane techniques was de- 
veloped. This analysis allows the PM stepping motor to be analyzed for failure to re- 
spond to given stepping rates with various values of load and friction. The nonlinear 
analysis was simplified by assuming that the stator windings are supplied by a constant 
current drive circuit that could compensate for the electrical (L/R) time constant. 

Thus, the results represent an upper bound on stepping motor performance. 

Based on the nonlinear analysis, it was found that the PM stepping motor cannot re- 
spond to a step command when the applied load torque becomes greater than 0. 707 of the 
motor’s stall torque. The PM stepping motor cannot follow a sequential set of step com- 
mands if, during the sequence, the rotor lags the command position by more than 
two steps. 

This result was extended to develop a set of dimensionless curves which express 
maximum normalized stepping rate as a function of normalized damping and normalized 
applied load torque. These curves allow the maximum stepping rate of a permanent mag- 
net stepping motor to be determined from simple measurements of its parameters or 
from manufacturer’s data. 
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Because the dimensionless curves assume a constant current source, the maximum 
stepping rate will be decreased if the motor is driven from a voltage source. The step- 
ping rate in this case will be decreased as a result of the stator winding electrical (L/R) 
time constant and induced electromotive force, both of which cause the rotor position to 
further lag behind the intended command position. 

Stepping motor response to a single-step command is usually highly underdamped. 
However, the maximum stepping rate curves show that improving the single- step re- 
sponse by adding damping will lower the maximum stepping rate that can be achieved. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 19, 1968, 

122-29-03-01-22. 
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APPENDIX A 


SYMBOLS 


B 

magnetic flux density, T 

n rt 

constant relating rotor position 

D 

viscous damping of motor and 
2 

load, (g)(cm )/sec 

in degrees 6 to mechanical 
degrees of rotation 

D 

normalized damping 

n 

number of solenoid turns 

E 

constant value of stator input 

R 

stator resistance per phase, 


voltage, V 

S 

Laplace operator, 1/sec 

E(S) 

Laplace transform of p(t), V 

t d 

developed motor torque, N-m 

E(t) 

stator input voltage, V 

T f (t) 

frictional torque, N-m 

E v 

stator induced electromotive 

Tf 

normalized frictional torque 


force (emf), V 

T L (f) 

applied load torque, N-m 

I 

constant value of stator input 

t l 

normalized applied load torque 


current, A 

T 

maximum motor torque per 

I(S) 

Laplace transform of l (t) 

max 

phase, N-m 

I(t) 

stator input current, A 

T s 

motor stall torque, both phases 

J 

inertia of motor and load, 

energized, N-m 


(g)(cm 2 ) 

t 

time, sec 

K I 

2 

solenoid constant, (£2)(sec)/m 

V 

normalized velocity 

K m 

bar magnet constant, m/(£2)(sec) 

X(t) 

variable 

k t 

motor torque constant, (N-m)/A 

a 

angular measurement variable, 

K v 

constant relating induced voltage 


deg 


and stator flux density, 

P 

line segment, cm 


(V)(see)/deg 

A() 

change in a variable about an 

L 

stator inductance per phase, H 


operating point 

if[] 

Laplace transformation 

6 

center of circular arc, deg 

l 

length of bar magnet, cm 

? 

damping ratio 

M 

magnetic moment, (A)(m ) 

e(s) 

Laplace transform of 0(t) , spe- 

m 

bar magnet pole strength, Wb 


cified in degrees except where 


noted in radians 
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0(t) rotor position defined so that 
1/4 rotor tooth pitch equals 
90°, specified in degrees ex- 
cept where noted in radians 

# c rotor command position, spe- 

cified in degrees except 
where noted in radians 

t normalized time 

At change in normalized time 

l/Ar normalized stepping rate 


<p(t) rotor position in mechanical 
degrees 

natural frequency, rad/sec 
(0) initial condition of a variable 

Subscripts: 

0 value at operating point 

1 stator phase 1 

2 stator phase 2 
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APPENDIX B 


LAPLACE TRANSFORM DERIVATIONS OF LINEARIZED FUNCTIONS 


The Laplace transformation can be applied by considering the following development. 
For a variable X(t) about an operating point, 


X(t) 

= X Q + AX(t) 

(Bl) 

dX(t) _ , 

+ Tdx(tj| 

(B2) 

dt 

- 1 
'S 
1 

T3 


d 2 X(t) _ 

/d 2 x\ , A d 2 X(t) 

(B3) 

dt 2 

\dt 2 / 0 dt 2 _ 



and so forth. Taking the derivative of equation (Bl) and assuming that Xq is a constant, 

— X(t) = iL (x n ) + A AX(t) (B4) 

dt dt u dt 

Combining equations (B2) and (B4), 



Taking the Laplace transform of equation (B5), 





(~) 

= && [AX(t)] - AX(0) - dt - 


(B6) 


In the same manner, taking the derivative of equation (B2), 



Combining equations (B3) and (B7) results in 





if[AE 2 (t)] - AE 2 (S) (B13) 

■S?[AI 2 (t)] = AI 2 (S) (B14) 


58 



1 die 


'dir 


r 

) A 

”dl 2 (t) 

1 


dt 

1 


— 

J 


.= S AI 2 (S) - AI 2 (0) 


dt / v»/ ft 

2 = s ai 9 (s) 2 

o ^ C 


-2'[AI 1 (t)] = AIj(S) 




dl^t) 


dt 


V dt 


» = S AI^S) - AI 1 (0) 


= S AI 1 (S) 


The derivative (dI 2 /dt) Q in equation (B15) can be determined by noting that since 
(d0/dt) Q = 0, the induced voltage E y Q = 0. The voltage supplied to winding 2 beet. 


AE 2 (t) = AI 2 (t)R + A 


dl 2 (t) 

dt 


L 


ae 2 (s) = (R + ls)ai 2 (s) 

Assuming that the stator input voltage to winding 2 is a step function results in 


or 


I 2 (S) = 


E 2,0/ 1 \ 

S \R + LS / 




dlo(t) 


dt 


= si 2 (s) = E2?0 

£ R + LS 


When the initial value theorem (ref. 14) is applied, 

limit ^2^ _ limit SE 2,0 = E 2,0 _ ^2^ 

t-0 dt R + LS L \ dt / 

0 


(B15) 

(B16) 

(B17) 

nes 

(B 18 ) 

(B19) 

(B20) 

(B21) 

(B22) 
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Thus, the Laplace transforms of the equations for the excursions about the operating 
point at the beginning of the step become 



T T 0 

+ — ^ AIjCS) - AT l (S) (B23) 

x i,o 


E. 


ae 2 (s) + ^0 _ 


AI 2 (S) = 


K 


RT 


- S A6(S) 


TC 2 T 2 
_T“1, 0 


R + LS 


(B24) 


AljOB) = 


— k’ - P S A0(S) 

N RT K T 1 1, 0 

R + LS 


(B25) 


Combining equations (B22) to (B24) gives the resultant linearized transfer function about 
an operating point at the beginning of a step: 


A0(S)=- 



*L,0 

1C 2 T 2 
^ 1,0 *- 


E, 


AE„(S) +— ^ 
1 S 


- (R + LS)AT l (S) 


JLS"4DL + JR) S 2 4^R +%5j + LKTll Q , 


N 


RT 


N- 


RT 


I n rt n rt 


Li5 - 


K 2 T 2 

“Tri,o 



(B26) 
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APPENDIX C 


DIMENSIONLESS PARAMETER CURVE COMPUTER PROGRAM 

A digital computer program was written to develop the dimensionless curves using 
the multistep phase plane analysis. The maximum stepping rate was determined by step- 
ping the motor a number of steps at a constant stepping rate. The phase trajectory was 
calculated and checked to see if it followed the command without missing a step. The 
stepping rate was increased until the motor failed to follow a command. 

The computer program developed is shown in table I. The program was developed by 
using the delta method for computing the phase trajectories. The delta method allowed 
the trajectory to be calculated by using small incremental sections. This method was 
also convenient for calculating time along the trajectory. The computer input variables 
were 

(1) Initial velocity, VTN 

(2) Initial position, THIN 

(3) Initial value of 6, DELTIN 

(4) Total time for the number of steps, TTOT 

(5) Total time between steps, DTAUS 

(6) Position increment used in computing trajectory, DTH 

(7) Normalized motor damping parameter, D 

(8) Normalized load torque, TL 

Twenty steps were commanded for each stepping rate. These were enough steps to 
determine whether the motor could follow the stepping rate with the given D and T t . 
Motor failure was determined by checking for the rotor velocity to change sign (motor 
reverses its direction) before the desired step command was reached. The stepping 
rates were varied in increments of 0.01, and the values of position were compared along 
the trajectory. The normalized time was computed, and a new value for delta was com- 
puted for each time the position of the rotor changed by a predetermined increment. 

This increment affects the overall accuracy of the trajectory and of the calculated time. 
Thus, it was made an input variable (DTH), and its effect on the response was checked. 

A value of DTH = 0.05 was found to give results with sufficient accuracy for the program. 

The maximum stepping rate was computed, using this program, for values of 
D = 0. 25, 0.5, 1.0, 1.5, 2.0, and 4. 0 and for values of T L =0, 0.2, 0.4, and 0.6. 
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TABLE I. - DIGITAL COMPUTER PROGRAM 


100 FORMAT (8E1 0.5) 

101 FORMAT ( 1 H 1 , 5X , 3H V I N , 9X , 4HTH I N , 7X , 6h DELT I N , 9X , 4HTT0T , 8X , 5H DTAUS , 
10X1,3HDTH,1 IX, 1HD, 1 IX, 2HTL/(8E 1 3 . 5) ) 

102 FORMAT (3E20. 6) 

103 F0RMAT(1H0,5X,8HVEL0CITY, 14X , 5HTHETA , 1 5X, 3HTAU) 

500 READ (5,100) VI N,TH I N,DELTIN,TTOT, DTAUS, DTH,D,TL 

WRITE (6, 101)VIN,THIN,DELTI N,TTOT, DTAUS , DTH, D, TL 

write( 6 , 103 ) 

TAU=0.0 
S I GN=1 .0 
V 1 =V I N 
TH 1 =TH I N 

WRITE (6, 1 02) VI ,TH 1 ,TAU 

DELTA=DELTI N 

N=1 

M=1 

1 0 RADS0= (V I*V 1 + (DELTA+TH 1 )* (DELTA+TH 1 ) ) 

TH2=TH 1+DTH 

TEM1=(DELTA+TH2)*(DELTA+TH2) 

I F (TEM1 . LE .RADSQ.) GO TO 20 
TH2=TH2-DTH 
DTH=- DTH 
S I GN=-S I GN 

TEM1 = (DELTA+TH2)*(DELTA+TH2) 

20 V2=S I GN*SQRT (RADSQ.-TEM 1 ) 

IF(TEM1 .GT.RADSQ) GO TO 21 
DTAU=2 . 0* (TH2-TH 1 ) / (V2+V 1 ) 

21 TAU=TAU+DTAU 
V 1 =V2 
TH1=TH2 

K=M-(M/5)*5 
I F(K)60, 60,70 
60 WRITE(6, 102)V1,TH1,TAU 
70 M=M+1 

if(tau.gt.ttct) go TO 500 

30 TEM=TAU - FLOAT (N— I )*DTAUS*4 . 0 



TEM2=TEM/DTAUS 





1 F(TEM2.GT .4.0) 

GO 

TO 

40 


GO TO 50 




ko 

N=N+1 





GO TO 30 




50 

IF(TEM2.GT.3-0) 

GO 

TO 

14 


1 F (TEM2.GT.2.0) 

GO 

TO 

13 


1 F (TEM2.GT. 1 .0) 

GO 

TO 

12 


11 VDV=COS (TH 1 )-D*V 1 -TL 
GO TO 15 

12 VDV=S I N (TH 1 )-D*V 1 -TL 
GO TO 15 

13 VDV=-COS (TH 1 )-D*Vl -TL 
GO TO 15 

14 VDV=-S I N (TH 1 )-D*V 1 -TL 

15 DELTA=-VDV-TH 1 
GO TO 10 

END 
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